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We have described electron spin dynamics in the presence of the spin-orbit interaction and disorder 
using the spin-density matrix method. Exact solution is obtained for an arbitrary 2D spin-orbit 
Hamiltonian and arbitrary smoothness of the disorder potential. Spin current depends explicitely 
on the disorder properties, namely the smoothness of the disorder potential, even in the ballistic 
limit when broadening by scattering is much smaller than the spin-orbit related splitting of the 
energy spectrum. In this sense universal intrinsic spin current does not exist. 
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Spin-orbit coupling brings about a number of interest- 
ing effects, one of which is generation of a spin flux in the 
plane perpendicular to the charge current direction. This 
phenomenon occurs in the paramagnetic system and has 
been very well known for quite a long time, see Refpi. It 
results from the fact that in the presense of spin-orbit 
coupling the scattering by impurities has an asymmet- 
ric character—. Extrinsic contribution exists only beyond 
the Born approximation in the scattering amplitude and 
leads to an accumulation of the spin density near the 
sample surfacei. 

It has been recently claimed^ that an analogous phe- 
nomenon can exist even without scattering by impurities, 
i.e. in the ballistic regime, the corresponding contribu- 
tion being called intrinsic. Later several papers appeared 
where the effect of scattering by impurities was taken into 
account for the Rashba mode&SiLSiSiiS, besides, the ex- 
trinsic contribution was describedii*i^ to explain the re- 
cent experiment done by Awschalom's groupie, see also 
experimenfeii. 

Here we study the dependence of the intrinsic contri- 
bution to the spin current on the form of the Hamilto- 
nian and properties of the disorder potential using the 
spin-density matrix methodii^. We have considered an 
arbitrary Hamiltonian, for example, the generalized 2D 
Rashba Hamiltonian with an arbitrary momentum de- 
pendence of the spin-orbit term or the 3D Luttinger 
Hamiltonian, and the arbitrary smoothness of the disor- 
der potential. In the case of the generalized 2D Rashba 
model and the arbitrary smoothness of the disorder po- 
tential an exact solution for spin current in the ballistic 
limit At 3> 1 is obtained, where A is the spin splitting of 
the electron spectrum and r is the transport scattering 
time. Spin current depends explicitely on the disorder 



properties, namely, the smoothness of the disorder po- 
tential, even in the ballistic limit. In this sense universal 
intrinsic spin current does not exist. The problem of 2D 
holes with the spin-orbit term has been numerically 
considered recently b y se veral groupsiSiii. Analytically 
it was studied in RefHa for the (5-scattering, where the 
incorrect answer was obtained due to the wrong calcula- 
tion method. The case of the arbitrary smoothness of the 
disorder potential has been analytically considered very 
recently in Refsii2i22i2i, some results were also reported 
in Refits. Below we will compare the results obtained 
here with those from Refl20l I also refer the reader to 
the introductory part of that paper where a complete list 
of references on spin current is presented. 

We consider the generalized 2D Rashba model which 
can be applied, for instance, for the case of 2D holes. The 
Hamiltonian of the problem is 
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H{p)^^ + ^a-nip), e,4p) = |- + Map^, (1) 

where a is the spin-orbit coupling constant, fl^ = 
cos(iV0), Qy = sin(iV(/)), flz — 0, eM{p) the eigenval- 
ues, M — ±1/2 the helicity values. The eigenfunctions 
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M=±l/2 A'=±l/2 



where D^^IP {Cl) is the rotation matrix^, 4> the angle of 
p, N the winding number, and the eigenfunction of 
the CTz operator. 

We will calculate the qy^ component of the spin cur- 
rent. This quantity is defined as 



qy, = aA^Tr J ^fip)^{S,yy + VyS.)^-l J _g_^(/^_(p) + /_^(p)). (2) 
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Here /(p) is the Wigner spin density matrix, Vy the y- component of the velocity operator and Sz — {\/2)(Jz the 
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spin operator. The last expression in Eq.|21) is given in 
the heHcity basis. The different components of the spin 
density matrix have the following relations to the average 
spin components: 



<S, >cx (/+-+/-+), 
<S-[n-n] >cx (/+_ - /_+), 

<s-n>^ (/++-,/—). 



(3) 



where n is the unit vector normal to the 2D plane (z- 
axis) . 

The general expression for the quantum kinetic equa- 
tion in the case of spin-orbit interaction, when the Hamil- 
tonian and the Wigner distribution function are matri- 
ces over the spin indexes, was derived, for example, in 
Refsii^*i^. In our case this equation is simple and reads 



dfjp) 

dt 



eE 



g/(p) 
9p 



-[^(p),/] = ^i{/(p)} (4) 



Here e is the charge of the carriers. The last term on the 
left hand side is the commutator and the expression for 
the collision term is given below. Now we write Eq.l^J 
in the hclicity basis where the Hamiltonian is diagonal. 
When doing that, we should take into account the fact 
that eigenfunctions X-fl/p depend on the direction of the 
momentum p, thus the matrix elements of the derivative 
df/dp in this basis do not coincide with the quantities 
dfMM'/dp 




df. 



MM' 



MM' 



dp 



'^[a, f]MM' 



^MM' 



dx 



Mp" 



M'p 



dp 



We see that there appears the commutator of the vector 
matrix a with /. Thus for Eq.Q in the linear response 
regime (E || x) we obtain 



eE cos ^^^^5 MM 



--iviivi '^-Y'~^^(-l'M'M'iP) - fAliliP)) + ji^Mij)) - ^M'{p))fMM'{p) = St{f{p))MM' 



(5) 
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Here is the equilibrium Fermi function corre- the helicity basis has the for: 

sponding to the helicity value M. The collision term in 
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Stif{p))MM'= I J2 {[^(eM,(pi)-eAf(p))+5(eM|(pi)-eM'(p))]i^jJ|/^i;;KpJ-/Ma/{(Pi) 



(27r?i) 



MiM'-^ 



-SieMAp) - ^M[iPl))[K'^qlqi^ppi) ■ /AfiM'(p) + JmmAp) ■ ^MJMj (^PPi)]}> 

I 



(6) 



where the kernel in the Born approximation in the scat- 
tering amplitude is: 



^m,m;KpJ = ^ic/(p-Pi)rn.i?iy,f^Kpji^5^^/^,^:Kp, 

(7) 

Here Ui is the 2D impurity density, J7(p — pi) is the 
Fourier component of the impurity potential. The quan- 



tities D 



'^Um (^PPi) depend only on the scattering angle 



D 



(1/2) 
1/2,1/2 



= D 



(1/2) 



D 



(1/2) 
1/2,-1/2 



1/2,-1/2 

n(i/2) 

-^-1/2,1/2 ^ 



= cos(iV6'/2), 
-ism{N9/2). 



(8) 



In this paper we consider only the case of small spin- 
orbit splitting, A/Ep ^ 1, A = app, and will be keep- 
ing only the terms linear in this small parameter. Note 



that expanding Eq.(|SJl with respect to this parameter, for 
the collision term in spin basis we obtain the expression 
which is identical to that given by Eq.23 of ref.l9^ We 
^ ) will be solving the problem in the helicity basis which is 
much more convenient and the solution can be obtained 
much easier. Moreover, further we will consider only the 
ballistic case, A ^ S/r, when one expects the intrinsic 
value for the spin current. In the case when the impu- 
rity potential has the azimuthal symmetry, it is possible 
to obtain a solution for the scattering potential with an 
arbitrary correlation length R. The most simple cases 
correspond to the following situations: 1) ma/Ti <C 
2) 1/R ^ ma/h, a — app'^ , Ti/ma is the spin-orbit 
length. We will see that these two cases correspond to dif- 
ferent answers for the spin current. The first case means 
a relatively short-ranged potential and was considered in 
Refs .l9.20> Note that it includes the limit of small angle 
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scattering when ma/Ti <C <C fc^- 

As here we consider the balhstic case, A ^ ?i/t, the 
problem can be greately simphfied, namely, in the col- 
lision term, Eq.®, we can neglect all the nondiagonal 
components of the spin density matrix since only diago- 
nal components can be proportional to r, which follows 
from Eq.|Sl). Moreover, since on the left hand side of 
Eq.(|3Jl there are only the first harmonics of the angle, 
the following solution obeys the system of Eas. H5l6() : 

f++{v) = f++{p) cos (t), /__(p) = /__(p)cos0, (9) 
/+-(P) = /+-(p)sin0, /_+(p) = /_+(p)sin<^. (10) 



It can be immediately seen from Eq.ljBjl with the help of 
the fact that the scattering kernel 



Ti ■ 

w{Q) = i^z\um\ uie) = u{\pi -p\) = 

Ui^Jp^ +pl~2ppi cose) (11) 



is the even function of the scattering angle 9. Indeed, 
from Eas. H5l6|l we obtain: 



eE 



(0) 



eE 



dp V+{p) 
gf(0) 



/++(?>) -T7^(«2/++(p)-«3/-(p+))-m25 j ^^sin2(^)(l-cos0)/++(p)U.o(12) 



V-{P+) 



dp V- {p) 



/__(p) + -P—{a^U+{p^) - a^f^^{p)) + m^a 



V+{P-) 



^sm^ — 1 - COS0 /++ p U=o 13 

ZTT dp 2 



.iV^(/0 
p 



-J++{p-)- 



6_)(/+_(p) -!-/_+ 
P+ 



ip))^ 
f — (P+)] - ^"^^a 



-ia4 



V+ip) 

de dw{e) 

2ti dp 



f++{p) 



:f-{p) + 



v-ip)- 

sm9sm{N9)f++{p)\ 

a=0- 



(14) 



In writing these equations we used Eqs.(|HIl- The coeffi- 
cients Ui in the above equations are given by the expres- 
sions 



de 

ai= I — VF(6') (cos 61- 1)(1 -I- cos 7V6I), 



"2 = I ^W{e){l-coaNe), 
J In 



de 

2^ 



1^(61) cos 6'(1 - COS iVe), 



de . 



04 = / — 14^(6*) sin 61 sin iV6', 

ZTT 



«5 = f ^wie){i-cos 

J ZTT 



9). (15) 



Note that here the scattering angle has an arbitrary value 
and is not assumed to be small. /+(p),/°(p) are the 
equilibrium Fermi functions which correspond to the he- 
licity ±, V±{p) — p/m ± Na/2 are the velocity values 
for a given p for ± bands. The other quantities entering 



Eas. p2ll4|) are defined in the following way 



p±=p±ma, V-{p-f-) 



P_ 

m 



:(iV-2), 



V+{p^)^^ + ^{N-2), & 
m 2 



app \ (16) 



The quantities y_(p+), Vj^[p-) are the velocities in the 
bands — , -I- for the momenta p+, p_, see Fig.l. These 
values of the momenta are connected in Eas. (|12ll4|l be- 
cause of the energy conservation under the elastic scat- 
tering. In writing Eas. (|16|l we took into account that 
ma pf, where p_F is the Fermi momentum. 

Eas. (|12ll4|l are written for a given value oi p. We see 
that the elastic scattering by impurities leads to admix- 
ture of the spin density matrix components which cor- 
respond to p+, p- values of the momentum. Different 
combinations of the ratio of the momentum and the ve- 
locity values likep/V+(p) orp+/VL(p+) are just different 
density of state values corresponding to different points 
in the momentum space involved in the transitions, see 
Fig.l. From the form of the coefficients in Ea. l|15|) we see 
that different terms in Eas. (|12ll4|l have a simple physical 
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FIG. 1: Schematics of the ± energy bands. Momenta p,p±,p]^ 
are shown, see the text. 



meaning, namely, they are just "transport" kinetic coefii- 
cients multipled (weighted) by the corresponding overlap 
factors sin2(iV6'/2), cos^{N9/2) or sin(7V6'/2) cos(iV6'/2) 
depending on the type of transition, i.e. whether it is 
intraband or interband scattering, see Eq.® and Fig.l. 
The very last term on the right hand side in each of the 
equations H12I14|I has a somewhat different origin. We 
can see that the scattering transitions within the band 
(+ or — ) and between the bands correspond to different 
values of the pi momentum entering the kernel of the 
collision term, see Ea. (|ll|) . In the case of the intraband 
transitions pi = p, while for the interband transitions its 
value is either or Thus, in the linear order with 
respect to a we obtain: 



W{p^p±)-W{q) = ± 



dW{q) ma 
~~dp 2~' 



W{q) 



W^(2p|sin-|). 
(17) 

Note that the quantity W{q) is the value of the scatter- 
ing kernel for the case of intraband transition when the 
modulus of the momentum is conserved, pi — p, and the 
modulus of the transferred momentum is q = 2p\ sin 0/2 1. 
Hence, we have taken into account the correction to W{9) 
in the corresponding terms of Eas. 1)12114(1 . keeping only 
the correction of the first order in a, i.e. all the pre- 
fectors in the terms which already contain the derivative 



dW{q)/dp are taken at a = 0. See, for example, the 
last term in Ea. (|12|l . where we also used that at a = 
f++{p) = f--{p) and given by Ea. (|^ . see below. 

All the quantities entering Eas. 1(12114(1 can be ex- 
pressed through the corresponding values at p. Writing 
Eas. l(12ll3() for the momenta values p±, and using the 
relations flip) = f_{p+), fl{p^) = f_{p), we obtain 



I—{P+) = f++{p) 



/++(p_) = /__(p) 



maNC 
Ap 



maNC 



Ap 



+ 



™«/++(p)|a=o dai _ 
2A dp 

maf++{p)\a=o dai , 



2A dp 

+||(iv-ua8) 



where E is given by Eg. 1(21(1 and A,C are expressed 
through coefficients a; by the formulas: A — (ai — a2 — 
03)72, C = ai + 05 



/df) 
— W (e) {cos 9 COS N9 - I), 
2n 

[ ^Wi9){cos9 -1) cos N9. 
27r 



(19) 



While deriving Eas. lfTHji . we used the expansions 

Wip- 



■P-] 



W{q)-ma^^, 
dp 

W{p+ ^p+) = W{q) +ma^^^^, 



(20) 



where the right hand side of these equations is written 
for the momentum value p. Note the difference in the 
coefficients in the Eas. ipU)! compared to Ea. ((T7jl . Using 
equation 



[/++(P) + /-(P)]U=0 



.dfo 



-^,E^eE^, (21) 
ma5 op 



which can be easily found from Eas. ((12113(1 . with /o being 
the equilibrium Fermi function at a — 0, we find the 

final expression for the quantity in question, (/,( {p) + 

f \-{p)), in the case ma/fi <^l/R: 



P 2a5 



04 fdai 2N 2(N- 1) 
A \ dp p p 



da/^ 
dp 



(e+-e-) 



iU-ip) + f-+{p)) 



(22) 



r 



where Ui are defined by Ea. p5(l . Let us check some lim- 1) We start with the simplest case of the (5-correlated 

iting cases from Ea. (|^ . potential, when scattering is isotropic, i.e. W{9) is a 
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constant. Then at iV ^ 1 the coefhcient 04 = 0, and 
from Ea. (|22ll we immediately obtain a simple result: 



p p il 

(23) 

With the help of Eas. H2llU|) for the spin conductivity at 
7^ 1 we obtain: 



eN 
"8^ 



(24) 



Note that this result also follows trivially from the di- 
agrammatic calculations since the vertex correction is 
identically zero. It happens because of a very simple 
reason, namely, the vector vertex contains the first har- 
monics but the Green's functions only the third one (at 
N — 3), and their overlap is zero, as it happens here with 
the coefficient a^. 

2) N=l. We know^iSiLSiSiiS that in this case for an ar- 
bitrary scattering (not nesessarely small angle or isotropic 
one) we should obtain a zero value for qy^. Indeed, in 
this case we have 04 = —A, ai = —04, 02 — 05 and from 
Ea. (|22|l we easily obtain: 



The result we obtain, i.e. Eg . H28jL coincides exactly with 
the corresponding result of Ref/23. 

4) Finally, we consider the case of a very smooth scat- 
tering potential when the impurity radius obeys the in- 
equality 1/R <C ma/fi ^ kp.'^^ Physically it means 
that elastic transitions between different bands (with op- 
posite helicity) are forbidden since the quantity W{9), 
see Ea. (|ll|) . is exponentially small for these transitions. 
Then in Eas. (|12ll4|l we should drop all the terms contain- 
ing p± quantities (and also derivatives dW/dp), and with 
the use of the fact, that in the case considered <C 1, the 
final result reads: 



p op op 



(/+_(p)+/_+(p)). 



(29) 



Integrating in Eq.Q between the points pp ~ - 
Pf, see Fig.l, we obtain for the spin conductivity 



eE 
P 



eiV 



(30) 



P 



dp 



if+-ip)+f-+{p))- 
(25) 

Expanding the first function on the left hand side of this 
equation with respect to a, we obtain a zero value for the 
spin current. 

3) Now consider the small-angle scattering case, <C 1. 
Then from Eq.^ it follows : 



NeE 
P 



(iV2 _,. 1) as dp 



rndtE 2N{N^ - N + 1) 
p (iV2TT) 



{f+-{p) + f-+{p)m 



Finally we need to calculate the quantity {l/ac,){das/dp) 
in the limit of small-angle scattering. For this quantity 
we obtain 



1 da.t; 



C-oo dq dW(q) 2 
^ _ 1 Jo 27r dq ' 'i' 'i 

^^dp - p /-|£M/(q)q2 



3 
p 



(27) 



where q = p\9\. With the help of Eq.(j2Z|), we finally 
obtain for the spin conductivity 



eN . {-5 + 2N -2N^ + 3N^). 

^ 's^^I' + ^v^TT^ 1 = 

eN {N"^ ~1){N -2) 



(iV2 + 1) 



(28) 



The first term (unity) inside the square brackets origi- 
nated from the first term on the left hand side of Eq. (|26|l . 



which for A'^ = 3 is twice as large compared to the S- 
scattering casei^^ 

In conclusion, we have studied how intrinsic spin cur- 
rent depends on the form of the Hamiltonian and the 
scattering potential properties. As an example, we have 
investigated the generalized 2D Rashba model when the 
spin-orbit term contains an arbitrary dependence on the 
electron momentum. We have found an exact analytic so- 
lution for the intrinsic spin current in the ballistic limit, 
when spin-orbit splitting is much larger than the disor- 
der induced broadening. In contrast to the case of a 
common, linear in momentum Hamiltonian, the intrin- 
sic current does not necessarily vanish. However, even 
in the ballistic limit indicated above, its value depends 
explicitely on the disorder properties. More precisely, 
the result, being independent on the spin-orbit coupling 
constant, is different for a different correlation radius of 
the disorder potential. The situation is similar for the 
other Hamiltonians. For example, we have also studied 
this problem for the Hamiltonian of 3D holesi^ In the 
diagrammatic language it means that even in the ballis- 
tic limit the vertex correction is always as important as 
the contribution of the empty loop and physically it is 
not correct to consider them separately. In this sense, 
universal intrinsic spin current does not exist. 
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